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Abstract 

We study the process B~ — > K~Xco considering intermediate charmed me- 
son rescattering effects. For this decay mode the naive factorization ansatz 
would predict a vanishing amplitude. We estimate contributions from the 
_£)(*)£)(*) _ ► K~Xco rescattering amplitudes, and compare the result with re- 
cent experimental measurements. We find that rescattering effects are able to 
produce a large branching ratio consistent with measurements by Belle Col- 
laboration. We also consider rescattering effects in B~ — > K~J/ip, arguing 
that they play a similar role in producing a large branching fraction for this 
colour-suppressed decay mode. 



1 Introduction 



Understanding strong interaction effects in weak exclusive heavy hadron decays is of great 
importance to gain information on fundamental aspects of strong and weak interaction 
phenomenology, both in the Standard Model and beyond. In this respect the factorization 
ansatz, that allows a treatment of nonleptonic decay amplitudes by factorizing hadronic 
matrix elements of four-quark operators as products of two current matrix elements, has 
been a widely used working tool for the analysis of B decays to charmed and charmless 
hadrons in the final state. In those decays in which the effective Wilson coefficients are 
not colour suppressed and the tree-level (V — A) x (V — A) current matrix elements do 
not vanish, it is found that factorization provides a reasonable description of data Q. 
However, the recent observation of the decay mode B~ — > K~Xco, reported by the Belle 
Collaboration together with the measurement of the branching fraction ||: 

B(B~ - K- Xc0 ) = (6.0t?i ± 1.1) x lO" 4 , (1) 

demonstrates the inadeguacy of the factorization model in the calculation of nonleptonic 
B decay amplitudes for colour suppressed B to charmonium transitions. A large non- 
factorizable term is needed to account for the observed branching ratio. As a matter of 
fact, the result ([!]) implies that the rate of B decays into a kaon and the ++ state of the 
charmonium system, Xco, is comparable with the B decay rate into a kaon and J/ift, and 
indeed the measurement of the ratio of the two branching fractions, reported by the same 
Collaboration, is: 

l[T-Z T-m = {om+ -™ ± °" 5 ± om) ■ (2) 

The experimental results ([J) and (0) are in conflict with the vanishing of the amplitude 
of B — > Kxco computed by the factorization ansatz, while the amplitude governing B — ► 
KJ/i/j is different from zero in the same approximation. This can be easily shown: the 
effective Hamiltonian governing both the transitions 1 : 

H w = % { v cb v: s (cMOM + cMOM) 

- (3) 

i 

involves only vector and axial-vector cc operators: 

0i = (cb) V -A{sc)v-A 

^^We neglect the B~ annihilation transition, which is governed by the CKM matrix element V u b- 
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2 = (sb) v _ A (cc) v _ A 

03(5) = (sb) V -A^2(qq)v-A\y+A] 
q 

^4(6) = [Sib^V-A^iqjQiW-AW+A] ( 4 ) 

q 

3 

°m = ^b) v ^ A ^e q (qq) v+A[v _ A] 
3 

^8(10) = ^{sibj)v-A^2e q {qjqi)v+A[v-A] 
1 q 

are color indices and (qq)v^A = <?7 M (1 T75)?)> an d therefore the factorized amplitude 
^ F (IT ^- Xc o) = %K b K;h+ £ a t ](K-\(sb) v _ A \B-)( Xd0 \{cc) VTA \0) (5) 

i=3,5,7,9 

vanishes since the current matrix elements (Xco|(3c)v,a|0) are zero. Instead, the experi- 
mental result (D corresponds to A exp (B~ -> K~x c0 ) = (3.39 ±0.68) x 1(T 7 GeV. On the 
other hand in the case of B — > KJ/ip, since (J/ip\(cc)v\0) ^ 0, a nonvanishing factorized 
amplitude, analogous to @, can be obtained once the combinations of Wilson coefficients 
0-2 = C2 + Ci/N c and a, = q + Ci + i/A^ c ,0 and the matrix element (K~\(sb)v-A\B~) are 
provided. 

Corrections to naive factorization involve gluon exchanges between the charmonium 
system and the quarks in B and K mesons. For a class of nonleptonic B — > M\Mi decays 
it has been argued that, in the large m& limit, non factorizable corrections are dominated 
by hard (perturbatively calculable) gluon exchanges, while soft effects are confined to the 
(B, Mi) system, where M\ is the meson picking up the spectator quark in B decay. This is 
the case of several processes where M\ and M2 are light mesons. However, when the meson 
which does not pick up the spectator quark is heavy, such a result no longer holds 0. In 
order to apply the QCD-improved factorization model to B decays to charmonium plus 
a kaon, either the cc state should be considered light with respect to the B meson, or one 
has to invoke the small transverse size of the cc system in order to assume a tiny overlap 
of the quarkonium wave function with the kaon wave function. However, an analysis 
of non factorizable corrections due to hard gluon exchanges in B — > Kxco has revealed 
the presence of infrared singularities, showing a difficulty of the method when applied to 
this decay mode ||. In the present note we investigate another effect, namely the K\ c q 
production by rescattering of open charm mesons Z) (*)£)(*) etc. primarily produced in B~ 
decays. The corresponding amplitude is mainly obtained by the operators 0\ and O2 in 
2 The definition of 0,7 and ag includes a factor e. 
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(|3|), and therefore this process could produce a sizeable contribution to B — > Kxco owing 
to the relatively large values of the corresponding Wilson coefficients C\ and c<i- Analogous 
effects were investigated in || , and have recently received new attention || [7| . 

Rescattering of intermediate D^D^ mesons can also contribute to the transition 
B — > KJ/ip and therefore we also consider this decay mode. Our conclusion is that, 
although the calculation presents uncertainties the size of which we shall try to assess, 
rescattering effects represent a non-negligible contribution to the decay channels B~ — > 
K~Xco an d B~ — > K~J/i/j. 

2 Process B~ -> D^~D^° K~XdO 

In the charm sector, rescattering effects have been recognized as a source of sizeable 
contributions in hadronic D meson decays. As the mass of the decaying B meson is 
larger than the D meson mass, one could suppose a minor role of such processes in B 
transitions, since one naively expects that high momentum final state particles move fast 
away from the interaction region without having the possibility to rescatter ||. However, 
in a number of analyses it has been shown that rescattering effects can play an important 
role even in B decays [§, [L(], [□], [T2| . 



It is worth attempting an estimate of the size of rescattering effects in color- suppressed 
B decays to final states containing heavy particles. We concentrate on two-body charmed 
meson contributions which can be included through a number of dynamical assumptions. 
We consider a set of amplitudes corresponding to the diagrams in fig.|I|, that represent 
t-channel contributions to the final state interaction. The charmed intermediate states 
and rescatter to Xco and K by the exchange of one resonance states, D and D*. 
We treat the exchanged resonances as virtual particles, with their propagators taken as 
Breit-Wigner forms. 

The analysis of the diagrams in fig.[3] involves the weak matrix elements governing 
the transitions B~ — > D^*'~D^*'°, and the strong couplings between the charmed states 
£)(*)-£)(*)o anc [ th e kaon an d Xco- There is experimental evidence that the calculation 
of the amplitude by factorization reproduces the main features of the B~ -> D^~D^° 
decay modes jnj. Therefore, neglecting the contribution of the operators O3-10 in (§), 



we can write: 



( D (.*)- D i*)o\ Hw \ B -) = 9z VcbVc * sai{D (*)°\( V _ Ar\B-}(D^-\(V - A),\0) (6) 



3 The role of inelastic effects in B decays has been emphasized in mI 
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Figure 1: Diagram contributing to the decay B — > K Xco- The boxes represent weak 
vertices, the dots strong couplings 



with ax — C\ + C2/N c . Using the heavy quark effective theory, the above matrix elements 
can be expressed in terms of a single form factor, the Isgur-Wise function £, and a single 
leptonic constant F |TJ]. This can be shown expressing the fields H a describing the 
negative parity J p = (0~, 1~) qQ meson doublet 



[PLY - Pal,] 



(7) 



in terms of operators P*^ and P a respectively annihilating the 1~ and 0~ mesons of 
four- velocity v (a = u, d, s is a light flavour index), and writing the B~ — > £)(*)° matrix 
elements as follows: 

< D°(v')\V^\B-(v) > = ^m B m D £(v ■ v')(v + v'f 
<D*\v',e)\Vi*\B-(v)> = -i^m B m D « £{v ■ v') e* p e^v^ (8) 
< D*°(v',e)\A^\B-(v) > = ^m B m D * £{v ■ v') e}[(l + v ■ v')g^ - v^v^} . 



In (U) e is the D* polarization vector 1 and ^(v ■ v') represents the Isgur-Wise form factor. 

The weak current for the transition from a heavy to a light quark Q — > q a , given at 
the quark level by g a 7 M (l — 75) Q, can be written in terms of a heavy meson and light 

iM 

pseudoscalars. The octet of the light pseudoscalar mesons is represented by £ = e / , with 

/ 



M 



2 



V 



\ 



7T 

K~ 



71 

-A/-7T 



K c 



K + 
K° 

1-2 



\ 



V J 



(9) 



and / ~ fn — 131 MeV, and the effective heavy-to-light current, written at the lowest 
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order in the light meson derivatives, reads: 

^ = ^[7^(1-75)^1] • (10) 

In this way, all the matrix elements < 0|g a 7 M (l — ^^)c\D^ [v) > are related to the single 
constant F: 

< 0\q a ^~f 5 c\D a (v) > = f Da m D y 

<0\q a Yc\D* a (v,e)> = f D *m D *e>* (11) 

with f Da = f D * = -J=. 

Other hadronic quantities appearing in the diagrams in fig.0 are the strong couplings 
£)(*)£)(*) /if and D^D^Xco- The D^D^K couplings, in the soft p K -> limit, can 
be related to a single effective constant g, as it turns out considering the effective QCD 
Lagrangian describing the strong interactions between the heavy D^D^ mesons and the 
octet of the light pseudoscalar mesons |IB|: 



£ x = % g Tr[H hl ^AlH a ] (12) 
with the operator A in given by 

A, ha = \{s}d^-id^) ha ■ (13) 

This allows to relate the Df'D^K couplings, defined through the matrix elements 

<D°(p)K-(q)\D*-(p + q,e))> = g DrD o K -(e-q) 
<D* (p, V )K-(q)\D:-(p + q,e))> = u*> a ep q^ g m - D . 0K - (14) 

to the effective coupling g: 

g 

9d*-d°k- = -2 v /m D m Z) * — 

JK 

g 

9d* b -d*°k- = 2y/m D .m D .— . (15) 

As for the coupling of the Xco state to a pair of D mesons, defined by the matrix 
element: 

(D°(p 1 )D (p 2 )\ Xc0 (p))=g DDXc0 , (16) 

an estimate can be obtained considering the D matrix element of the scalar cc current: 
(D(v')\cc\D(v)), assuming the dominance of the nearest resonance, i.e. the scalar cc state, 
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in the (v — t/) 2 -channel and using the normalization of the Isgur-Wise form factor at the 
zero-recoil point v — v'. This allows us to express gDD Xc o i n terms of the constant f Xc0 
that parameterizes the matrix element 

(0|cc|Xd)(?)> = f X co m XcO ■ ( 17 ) 

The method can also be applied to gD*D* X co- O ne obtains: 

m D m Xe0 _ n rn D *m Xc0 

9DD X co — ~ Z 7 9d*D* X cO — 1 7 • 

JXcO JXcO 

It is worth noticing, however, that the determinations of the couplings described above 
do not account for the off-shell effect of the exchanged D and D* particles, the virtuality 
of which can be large. As discussed in the literature, a method to account for such effect 
relies on the introduction of form factors: 

gi (t)=g i0 F t (t), (19) 

with g i0 the corresponding on- shell couplings (|14D, (|16|). A simple pole representation for 

2 ^2 



Af — m 



Fi(t) is: Fi(t) = 1 2 — ^_!_, consistent with QCD counting rules HI7| . The parameters 



in the form factors represent a source of uncertainty in our analysis. 

We have the elements for computing the diagrams in fig.|l]. The absorptive part of the 
amplitude (1) simply reads: 



^\{m 2 B ,m 2 Dp m 



327rm| 



B 



+ 1 



lmA 1 = -2 — j / dzA(B~ -> D*~ D )A(D*~ D -> K'Xco) (20) 



with A the triangular function. Analogous expressions correspond to the diagrams (2) 
and (3). Explicitely, the imaginary parts are given by 

Kf D *m D .y/m B m D (mB + m D ) . 1/2 , 22 2 x f 1 . /.x /.x 

ImAi = 2 327rm 2 m D ( m B^ m Dp m D> J ^z g DtDK {t)g DDx {t) 

_g0 _| 2 

( m 2 B -m 2 m + m 2 D \ m|,. 

\ 2m B m D J t — m 2 D 

Im^ 2 = 577 — 2 A ' {m B ,m D m Dt ) / dz g DaD * K (t)g D * D * x (t) 

SZTTmi, J-i 



' m 2 n — m 2 n + m 2 n * \ f \m 2 K — q ■ k 



5 ' 2m B m^ / I 



1 



1 + v ■ v D )q ■ k + v D ■ k(q° + v D ■ q)} 
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-[-(1 + v ■ v D )m D + v D ■ k(k° +v D ■ k 



lmA 3 = - K ^ Ds ? ^ r s m 2 mi?mD ( m B» m z) s * » m D* ) / ^dz g D * D , K {t)g D * D * x (t) 



B 

/ m| - m 2 ^ + m 2 ^ \ q% D ■ fc - k°v D ■ q 
* \ 2m B m* D ) t - m 2 D , ' ^ ' 



G Fyr rr * o m 2 B +m 2 K -m 2 A^ 2 ( 



m%, mi-, m 2% 



where: K = -^V^VLai, o = -, ol = — , together with: 

m B -k° , m B k°-m 2 Dgt m B q° - q ■ k m 2 B + m 2 D - m 2 D 
vv D = ,v D -k = >-, VD .q= } fc = — ! , 

m D . m Di m Di 2m B 

\ 1/2 (m 2 B ,m 2 D ,m 2 D ) 

\k\ = and m Da = m Da = m D *, m Ds 2 = m Ds ; m Dl = m D , m D2 = 

2m B 

m Dz = m D *. 

The dispersive parts of the amplitudes in fig.p] can be estimated using 

ReMO = l-PV [T ^^^T ( 22 ) 



7T 



« ■ s'-m 2 B 



with the thresholds sf^ given by: s$ = (m D * + m D ) 2 , sf^ = (m Ds + m D *) 2 and = 
(mo* + rriD*) 2 , respectively. It can be assumed that such expressions are dominated by 
the region close to the pole m 2 B . Therefore, we compute the integrals, that in general 
depend on the asympotic behavior of the spectral functions ImAi(s'), by using a cutoff 
not far from the B meson mass, chosen in the range 35 — 40 GeV 2 . 

A comment on other contributions to K~Xco via final state rescattering is in order, 
since also the D^*> and excitations could be considered as intermediate states. These 
terms are suppressed by smaller values of the universal form factors and of the leptonic 
decay constants; therefore, the amplitudes in fig.|l] represent the main contributions that 
need to be analyzed. 



3 Mode B~ -> K~ J '/ip 

Before attempting a numerical estimate of the amplitudes in fig.|l|, let us consider the 
decay mode B~ — > K~ J/ip. In this case the amplitude obtained in the naive factorization 
approach, keeping only the contribution of the operators 0\ and O2 in (H), is given by 

A F (B-^K-J/^) = 2^V cb V^a 2 f J/f mj / ^ K ( m 2 ^)(e* ■ q) 

= M^-Q) , (23) 
8 



with the constant fj/^ defined by 



(0\&fc\J/i>(j/ t e)) = fj/^mj/^ , (24) 

e the J/ip polarization vector, q the kaon momentum and one of the two form factors 
parameterizing the matrix element (K~ \ s^f^b\B~). 

Identifying eq.(P3"D with the experimental amplitude obtained from the measurement 
B(B~ -> K-J/ip) = (1.00 ± 0.10) x 10~ 3 [T§: A exp = (1.41 ± 0.07) x 10^ 7 , and using the 
value fj/$ = 405 ± 14 MeV, one obtains a result for the product \a 2 Ff K (m 2 j^)\. This 
determination of a 2 is mainly affected by the uncertainty on Fj 3 , scanning several form 
factor models, as done in one gets \a 2 \ = 0.2 — 0.4, while considering the calculation 
20|| one obtains |<Z2 1 = 0.38 ± 0.05. Such results are obtained using V c b = 0.040 and 



in 



V cs = 0.9735 that correspond to the central values reported by the Particle Data Group 



18 



As discussed at length in the literature, the above values of a 2 are different from the 
combination a 2 = c 2 + c\/N c of the Wilson coefficients in (Q). As a matter of fact, from the 
values d = 1.085(1.109) and c 2 = -0.198(-0.243) computed for m b (m b ) = 4.4 GeV and 
Ag)_ = 290 MeV in the naive dimensional regularization (or 't Hooft-Veltman) scheme 
glj, one would get: a 2 = 0.163(0.126). Q 

Therefore, nonfactorizable effects are sizeable in B~ — > K~J/if>, and indeed in a 
generalized factorization ansatz a 2 is treated as an effective parameter used to fit the 
data. The calculation in the framework of QCD factorization does not allow to reproduce 
the fitted value, although an improvement towards the experimental datum is obtained 
p2| . It is worth considering rescattering contributions of intermediate charm mesons, 
described by diagrams as depicted in fig.[| The hadronic information for determining 
such amplitudes are the same as in Section 2, with the only difference in the strong 
D^D^J/ijj couplings that can be expressed in terms of the parameter /jm, using the 



same vector meson dominance method applied to derive eq. ([18]). 



4 Numerical calculation and discussion 

In order to evaluate the amplitudes in figs.|l],^| we have to fix the values of the various 
hadronic parameters. The Wilson coefficient Oi, common to all the amplitudes, can 
be put to di = 1.0 as obtained by the analysis of exclusive B — > £)(*)£)(*) transitions. 

4 Similar values are obtained varying m^mj) and A-^-L. 
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(1) 



(2) 




(3) 



(4) 



Figure 2: Rescattering diagrams contributing to B — > K J /ip. 



240 MeV, in the range quoted by the Particle Data Group |I5 



Moreover, we use fo a 

and fn* = fn s consistently with our approach that exploits the large mq limit. For the 

Isgur-Wise universal form factor £, the expression £(?/) = ^ — * s com P a tible with the 

current results from the semileptonic B — ► decays. 

A discussion is needed about the D^D^*'K vertices. For the effective coupling g in 
(|i~5D one can use the CLEO result g = 0.59 ± O.OTi 0.07 obtained by the measurement of 
the D* width |23[ . This value is in the upper side of several theoretical calculations E3, |25 



We choose to be conservative, and vary this parameter in the range: 0.35 < g < 0.65 that 
encompasses the largest part of the predictions. 

We use the expression flT%| ) for the D^D^XcO vertices, with f Xc0 = 510 ±40 MeV ob- 
tained by a standard two-point QCD sum rule analysis. As for the couplings D^D^J/ip, 
expressions analogous to (|l8|) involve fj/^, for which we use the experimental measure- 
ment. To account for the off-shell effects of the exchanged particles, we use eq. (|l9D 
with two choices for the parameters: Aj = 2.5 GeV and Aj = 2.8 GeV, corresponding 
to typical values of the mass of the radial excitations of mesons. 

The results are reported in Table |l| One has to notice that the rescattering amplitudes 
contribute with different signs to the final result, with significant cancellations among the 
various terms. 

A few observations are in order. First, in the chosen range of values for the vari- 
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Table 1: Numerical results for the rescattering amplitudes 



B- - K~ X co 


Re„4 (GeV) 1mA (GeV) 


Ai (GeV) 




-(0.9- 1.7) x 10~ 7 -(0.5- 1.0) x 10" 7 
-(1.4-2.7) x 10- 7 -(0.6 - 1.2) x 10" 7 


2.5 
2.8 




Re.4 1mA 


A, (GeV) 




(0.1-0.2) x 10- 7 -(0.5-0.9) x 10~ 7 
(0.2-0.3) x 10~ 7 -(0.9- 1.7) x 10" 7 


2.5 
2.8 



ous parameters, the rescattering amplitudes are sizeable, and become comparable to the 
experimental ones. This observation can be made more quantitative. Assuming that 
the amplitude relative to B~ — > K^J/tp deviates from the factorized result because of 
the contribution of the rescattering term: A exp = A fact + A resc , one can constrain the 
values of Aj for the calculation of B(B~ — > K^Xco)- One obtains: A, ~ 2.7 GeV and 
B(B~ -> K- X co) = (1.1-3.5) x 10~ 4 , to be compared to fll|). The result seems noticeable, 
considering the rather schematic description of the rescattering process. 

The second observation is that a correspondence similar to that between B~ — > K~Xco 
and B~ — > K^J/ip is expected with analogous decay modes, namely B~ — > K~Xci 
and B~ — > K~Xc2, the mesons x&,2 being the axial vector and the tensor states of the 
charmonium system. It is worth observing that B~ — > K~Xc2 is another process the 
amplitude of which vanishes in the naive factorization model; in our approach, we expect 
a branching fraction analogous to that of B~ — >• K~Xco- 

We are aware of the various sources of theoretical uncertainty. The uncertainty related 
to the B~ — > £)(*)£)(*) vertices can be minimized gaining experimental information on 
these decay modes. The uncertainties in the strong vertices can be reduced by dedicated 
analyses using nonperturbative QCD methods (QCD sum rules or lattice QCD). However, 
all such uncertainties only affect the precise numerical predictions and not our main 
conclusion. We have found that rescattering amplitudes, describing a rearrangement of 
the quarks in the final state after the production of pairs of charmed mesons, not only 
cannot be neglected both in B~ — > K~XcQi both in B~ — > K~J/ip, but can provide a 
large part of the decay amplitudes. Analogous effects are expected to be important in 
similar colour suppressed decay modes, namely B~ — > K~Xci an d B~ — > K~Xc2- 
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